Comparison of Traditional Design 
Nonlinear Programming Optimization and 
Stochastic Methods for Structural Design 


Structural design generated by traditional method, optimization method and the 
stochastic design concept are compared. In the traditional method, the constraints are 
manipulated to obtain the design and weight is back calculated. In design optimization, 
the weight of a structure becomes the merit function with constraints imposed on failure 
modes and an optimization algorithm is used to generate the solution. Stochastic design 
concept accounts for uncertainties in loads, material properties, and other parameters 
and solution is obtained by solving a design optimization problem for a specified 
reliability. Acceptable solutions were produced by all the three methods. The variation in 
the weight calculated by the methods was modest. Some variation was noticed in designs 
calculated by the methods. The variation may be attributed to structural indeterminacy. It 
is prudent to develop design by all three methods prior to its fabrication. The traditional 
design method can be improved when the simplified sensitivities of the behavior 
constraint is used. Such sensitivity can reduce design calculations and may have a 
potential to unify the traditional and optimization methods. Weight versus reliability traced 
out an inverted-S-shaped graph. The center of the graph corresponded to mean valued 
design. A heavy design with weight approaching infinity could be produced for a near- 
zero rate of failure. Weight can be reduced to a small value for a most failure-prone 
design. Probabilistic modeling of load and material properties remained a challenge. 
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Objective 

• Summarize work performed supported by the Subsonic Fixed Wing project 
under NASA’s Fundamental Aeronautics Program. 

• Focus: 

• Comparison of Three Structural Design Methods 

- Traditional Design 

- Nonlinear Optimization Programming Techniques 

- Stochastic Design Optimization 

- Probabilistic Analysis - Four methods 
- Solution for a three bar truss 

• Stochastic Design Optimization Code (SDOC) 

- CometBoards testbed as Optimizer 

- FPI of NESSUS code as Probabilistic estimator 

- MSC/Nastran as Structural Analyzer 

• Results for Three Design Problems 

- Tapered Beam - Determinate 

- Clamped Beam - Indeterminate 

- Raked Wing Tip Structure of Boeing 767-ER airliner - Complex 

• Conclusions 
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Methods for Structural Design 


Traditional Design 

Find n design variables {x} within 
bounds to satisfy (fully utilize) 
m constraints gj <= 0 
Constraints are manipulated to obtain 
the design. 

A condition to minimize weight is not 
imposed 


Deterministic Optimization 

Find n design variables {x} 
within bounds 
to minimize weight W(x) 
subjected to m constraints g- <= 0 
Nonlinear mathematical programming 
algorithms 


Stochastic Design 

Find n mean values for n random 
variables {x} within bounds 
to minimize mean value of weight W(x) 
subjected to m constraint g- (p) <=0 HZ 
defined as: P{gj <=0} >= pj 
Same as deterministic design but the 
parameters are random in nature 


Probabilistic analysis 

Primitive random variables are defined 
through a distribution function with a 
mean value and a standard deviation 
Response is calculated via MPP of 
NESSUS code. The influence of each 
variable on the probability of constraint 
satisfaction can be obtained by 
calculating probabilistic sensitivities 
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Traditional Design 
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Basic problem: Calculate design variables d 1 and d 2 to find the weight 
subject to strength and displacement constraints. 

Traditional design method is based on a fully utilized concept 


Design is first obtained for strength limitation: 

f My\ ( My ^ 


<7„ = <7, 


and cr„ = cr„ 


V I J 


C/ 1 = > /6M 1 / ba o ; d 2 = ^6M 2 1 ba o 

Stiffness constraint can be written as: 

Pa 3 


9,(d v d 2 )- 


L+A 

C/ 3 c/3 
U 1 U 2 


Design can be adjusted for violated stiffness constraints: 
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Constraint is satisfied by many sets of (d 1 ? d2) an( l weight is not unique 
Weight is not used in traditional design but back calculated 
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Deterministic Optimization 


• Cast as a nonlinear mathematical programming problem 

• Find x that minimize W(x) subject to g(x) < 0 , h(x) = 0 and x lb < x < x ub 

where W is an objective, x is a vector of design variables, g is a vector of inequality 
constraints, h is a vector of equality constraints, and x lb and x ub are vectors of lower 
and upper bounds on the design variables. 

• Applications of nonlinear programming include: automobile design, naval architecture, 
electronics, computers, aerospace engineering, such as aircraft and spacecraft design, 
etc. 

• Implemented several methods in Stochastic Design Software Code (SDOC) 

- MSC/Nastran was the analyzer 

- Optimizers were taken from CometBoards 
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Organization of Stochastic Design Optimization Code 

(SDOC) at NASA GRC 




Can solve structures, airliner and jet engine design problems 
Being interfaced in OpenMDAO Framework at NASA GRC 
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Stochastic Design Optimization 

• Similar to deterministic optimization but parameters are random in nature 
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Stress x as a random variable 



• For stress allowable we require three assumed parameters: 

- Distribution function (Normal) for strength (random variable) 

- Mean value 

- Standard deviation 

• For stress response we assume a Normal distribution function 

• Probability of success also needs to be chosen 

• Accuracy of design depends on these parameters four of which were assumed 

Glenn Research Center at Lewis Field 



CD-08-83071 


Summary of Probabilistic Analyses Methods 


Perturbation method (PM) 

Direct Monte Carlo simulation (DMCS) 
Latin hypercube simulation (LHS) 

Fast probabilistic integration (FPI) 


Ten primitive random variables (qrqi 0 ) for 3 -bar truss: 
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Three bar truss 
Covariance matrix included bar areas with 
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Similar for other variables 
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Probabilistic Analyses Results 


Parameter 

Perturbation 

method 

(PM) 

Direct 

Monte Carlo 
simulation 
(DMCS) 

Latin hypercube 
simulation 
(LHS) 

Fast probability 
integrator 
(FPI) 


Mean 

value 

Standard 

deviation 

Mean 

value 

Standard 

deviation 

Mean 

value 

Standard 
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Mean 
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• Performance was satisfactory for all four methods 

• Monte Carlo simulation is numerically expensive while FPI took the least 

• FPI is used for probabilistic analysis throughout this study 
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Comparison of Results for Three 
Structural Design Methods 
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Design of a Tapered Beam 

Calculate depths d 1 and d 2 for strength and displacement limitations 


Variable 

Deterministic 

value 

Mean value 

Standard deviation 



Stress limit: <j 0 

20 ksi 

li = 25 ksi 

' CT 0 

o = 2.5 ksi 

Disp. limit: s 

1.0 in. 

p g = 1 .25 in 

o s = 0.125 in. 

Modulus: E 

10 ksi 

p E =10 ks i 

<7 E = 1 ksi 

Density: p 

0.1 lb/in. 

^=0.1 

o p = 0.005 in. 

Load: P 

20 kip 

p p =15 /c/p 

Op = 3 kip. 

Depth: di, d 2 

Calculate 

p d ^ = 25 in and p d2 =10 

<7 di = 3 in and a d2 = 2 in 

Limitation 


Assumed but not measured values 
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Design Solutions for Tapered Beam by Different Methods 


Method 

Design variables in 
inch 

Weight 

(Ibf) 

Max stress 

Max disp 

4 

d 2 



Proration Technique 

19.4 

13.8 

955.4 

11.4 

1.0 (active) 

Engineering Method 

25.0 

10.4 

1018.1 

20.0 

1 .0 (active) 

Optimality Criteria 

20.4 

12.5 

946.9 

13.8 

1 .0 (active) 

NLPQ Algorithm 

20.4 

12.6 

951.9 

13.6 

1.0 (active) 

Cascade Strategy 

20.5 

12.5 

951.9 

13.8 

1.0 (active) 

Stochastic Design for 1 failure 
in one million samples 

25.1 

10.2 

1028.7 

— 

active 
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Summary 


No unique solution even for the 
determinate beam 
Variation in solutions 

- Weight: 8 percent 

- Design variables: 35 percent 
Performance 

- Acceptable by all methods 
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Stochastic Design Solution for Tapered Beam 


N 

(One failure in N samples) 

Design variable di 
in inch. 

Variable d 2 

Weight in lbf 
(mean value) 

2 

17.2 

10.6 

801.7 

10 

19.8 

12.2 

921.0 

100 

21.5 

13.2 

1000.5 

1000 

22.7 

14.0 

1055.4 

10,000 

23.6 

14.6 

1100.6 

100,000 

24.5 

15.1 

1140.6 

200,000 

24.7 

15.2 

1152.1 

500,000 

25.1 

15.4 

1167.0 

1 million 

25.3 

15.6 

1178.0 

1.25 million 

25.4 

15.6 

1181.6 

2 million 

25.5 

15.7 

1189.0 
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Summary 

Weight increased when risk was reduced 
& vice-versa 

For 50% rate of success weight = 802 lbf 
Weight — ► to oo for zero risk 
Weight — ► to 0 for most unreliable design 
Extreme values can not be captured 
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Design of a Clamped Beam by Different Methods 


Calculate depths d 1 and d 2 for strength and displacement limitations 
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Variable 

Deterministic 

Mean value 

Standard deviation 


value 



Stress limit: a 0 

20 ksi 

p ao = 25 ksi 

=2.5 ksi 

a 0 

Disp. limit: 8 

0.25 in. 

p s = 0.251 in 

o 5 = 0.005 in. 

Modulus: E 

1 OK ksi 

p E =10 K ksi 

a E =1 K ksi 

Density: p 

0.1 lb/in. 

^ =0.1 

o p = 0.005 in. 

Load: P 

25, 50, 25 kip 

p p =25.01,49.98,25.01 

(t p = 0.01,0.01,0.01 kip 

Depth: di,d 2 

20.34,17.49 

p di = 23.805 in and 
p d2 =12.805 in 

cr d = 0.005 in anda d2 = 0.005 in 

Limitation 


Assumed but not measured values 
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Design Solutions for the Clamped Beam 


Method 

Design variables in 
inch 

Weight 

(Ibf) 

Max stress 
(ksi) 

Max disp in 
inch 

4 

d 2 



Proration Technique 

23.5 

12.9 

699.1 

15.13 

0.25 

Engineering Method 

21.5 

15.3 

707.1 

19.0 

0.25 

Optimality Criteria 

23.8 

12.8 

702.7 

15.4 

0.25 

NLPQ 

24.5 

13.5 

728.6 

14.2 

0.25 

Cascade Strategy 

24.5 

13.5 

728.6 

14.2 

0.25 

Stochastic Design for 1 failure 
in 1 0 samples 

29.5 

16.35 

881.9 

— 

— 
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Design of a Raked Wing Tip of Boeing 767-400 ER Airliner 
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• The structure shown was made of composite 
and different types of metals 

• Optimum design was obtained for 

1 failure in N samples, Max N = 2 Million 

•13 sets of design variables 

• 227 groups of strain & displacement 
constraints 

• Solution required 5 days of execution 
in a Linux workstation. 

• Limited information can be provided 
because data is Boeing proprietary 
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Ten Primary Components of Wing Tip 



"r 


Component 


68.5 


Tip rib skin with 
3 webs- Passive 
GRC Sub 7 
Boeing #29, 30-7 


Passive? 

5.31 Ibf (0.063/0.1) 

Pair of titanium doublers- T.l 
GRC Sub 8 
Boeing #12, 13 



The wing component (subproblem) 


Leading edge- 


Passive 



T.E. 

GRC Sub 3 
Boeing #28 


Center 

Passive | i GRC Sub 1 
Boeing #17 

LE. 

GRC Sub 2 
Boeing #27 
3 - Metal 

Root rib 


Active 





Frame structure 


Face sheets-composites 


- Composite upper and 
lower skin panel 
GRC Sub 9 and 10 
Boeing #10, 11 


9 Ibf (0.063) 
Leading edge skin 
GRC Sub 6 
Boeing #14 

Rear web- Active 
(metal plus composite) 
GRC Sub - 5 
Boeing #23, 24 


Front web- Active 
(metal plus composite) 
GRC Sub 4 
Boeing #15, 19, 16, 20 
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• Optimized four composite components (shown in red) 

• The six metallic components were passive (shown in blue titles) 
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Orthotropic Material Model 


Material properties 

Deterministic 

values 

Probabilistic values 

Mean value 

Standard deviation 

Material w Fabric 

8.1x10 6 psi 

8.505x1 0 6 psi 

0.638x10 s psi 

Young’s modulus, E 

Shear modulus, G 12 

7.1x10 5 psi 

7.35x1 0 5 psi 

0.551x10 5 psi 

Honey comb 

4500 

4725 

334 

Transverse shear 
modulus G 1z 

G 2z 

2500 

2625 

197 

Strain allowable 

4000 ps 

5000 

375 


Limitation: Mean values & variations were assumed not measured 
Failure theory: Maximum failure theory 
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Design Variable and Constraint Definition 


Design 

variable 

Location 

Honeycomb 
Thickness in 
inch 

Number 
of plies* 

1 

Section 1: Upper & lower 
skin panel 

0.65 - 0.75 

7 - 16 

2 

Section 2: Upper & lower 
panel 

0.65 - 0.75 

18-27 

3 

Section 3: Upper & lower 
panel 

0.65 - 0.75 

28 - 38 

4 

Section 4: Upper & lower 
panel 

0.5 

7 - 14 

5 

Section 5: Upper & lower 
panel 

0.5 

15 - 21 

6 

Section 6: Upper & lower 
panel 

None 

9 - 17 

7 

Section 7: Upper & lower 
panel 

None 

18 - 27 

8 

Section 8: Upper & lower 
panel 

None 

28-38 

9 

Web of both spars 

None 

8 - 15 

10 

Front & aux spars 

0.18-0.36 

!3 sq. in. 

11 

Rear spar 

0.27 - 0.343 

12 

Upper -passive 

36 plies 

13 

Lower -passive 

38 plies 


• Structure was made of (7-38) plies, honeycomb core and rod elements 

• Behavior constraints: 227 (203 strain for upper & lower panels + 16 for rods) 

• 3 translations (z-dir <= 13.68 in.; mean = 17.1, stdev = 1.71); 

• 1 rotation (5°; mean = 6.25, stdev = 0.625) 
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Design Variable Grouping for Panels 


Upper panel 




Lower panel 


DV 1 - It green 
DV 2 - blue 
DV 3 - green 
DV 4 - orange 
DV 5 - dk green 
DV 6 - red 
DV 7 - pink 
DV 8 - It pink 


• Design variable grouping was formulated by considering elements whose 
thickness range was similar. 
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Initial Analysis for Load Case Determination 


Nominal Response 

Load Case 

Max. 
Disp. in 
inch 

Micro-Strain in C- panel 

Strain 
Energy in 
in.-kip 

Lower 

Upper 

1 

-5.6 

1042 

3243 

9 

2 

11.7 

5744 

1713 

40 

3 

10.6 

5648 

1579 

34 

4 

10.4 

5689 

1512 

33 

5 

-5.5 

985 

2998 

9 

6 

13.7 

5897 

3773 

54 

7 

12.4 

5189 

3409 

44 

8 

-10.0 

1787 

4737 

28 


• Frequencies: 1 st mode =12.7 Hz; 2 nd =16.4; 3 rd =17.0 

• Lower composite panel carried more load 

• Weight of composite materials = 84.5 lbf 
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Deterministic Optimization Results 



Boeing Design 

Design Model 1 

Design Model 2 

Weight, lbf 

84 5 

70.966 

67 792 

Design variables ; 
in inch 




1 

0.0958 

00680 

0.0680 

2 

0.1791 

0.1871 

0.1 7 SB 

3 

0,2730 

0.3012 

0.3053 

4 

0.0899 

0 0680 

0.069-2 

5 

01516 

0 0680 

0.0600 

6 

00998 

0,1231 

0.0808 

7 

01857 

0.1563 

0.1754 

8 

02727 

0 1202 

0.1224 

9 

0.0963 

0.0786 

0.080S 

10 

03065 

0 2700 

0.2700 

11 

0.2438 

0.1800 

0.1300 

12 

0.3060 

0 3060 

0.3060 

13 

0.3230 

0.3230 

0.3230 

Active 

constraints 


Six strains = 4000 /is 

Si* strains = 4000 ps 

Displacement 

1 3,68 m 

13. 14 

13.35 

Rotation 

A 24" 

4 48 

4.72 

Frequency, Hz 

1636 

16-45 

20 33 


• Weight was lighter by 13.5 lbf or 16 percent of original design 

• Strain was about 4000 jlis for six of 13 groups of design variables 

• Frequency was 16.45 Hz against 16.36 Hz for initial design 
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Stochastic Design Optimization Results 


Probability 
of single 
constraint 
Failure 

N 

One failure 
in N samples 

Model-1 
Load-A 
Strength + 
stiffness 

Model-1 

Load-B 

Strength or strength 
and stiffness 

Model-2 

Load-A 

Strength and stiffm 

0.5 

2 

54.726 

52.658 

53.474 

0.1 

10 

57.054 

53.782 

54.453 

0.01 

100 

59.620 

55.429 

56.661 

0.001 

1,000 

62.398 

56.664 

58.861 

IE-4 

10,000 

64.844 

57.727 

61.322 

IE-5 

100,000 

67.076 

58.969 

64.512 

IE-6 

1 million 

78.626 

60.553 

77.170 

8E-7 

1,250,000 

80.272 

60.716 

79.028 

5E-7 

2,000,000 

84.607 

61.057 

83.427 


• Weight varies depending on the failure rate - the less likely the failure the higher 
the weight 

• Each row represents a different probability of single constraint failure 
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Weight (Ibf) 


Inverted S-graph - for Model-2 and Load A 
x-axis in log scale, weight in y-axis 



• Weight can be approximated by two linear segments 

• First segment is for active strain constraints only 

• Both strain and displacement are active for the second segment 
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Strain Distribution in Upper Panel 


Boeing design 
Max Strain = 5897 ps 
Total Weight = 84.5 lbf 

Deterministic design 
Max Strain = 4000 ps 
Total Weight = 70.97 lbf 

Stochastic design for one failure 
in 2 Million samples: 

Model- 1, Load A 
Max Strain = 4670 ps 
Total Weight = 70.354 lbf 


ived Results. Subcam? 3. Strain Tenaor^SCAPei (N ON-LAYERED) 4.000-03^B 

- ■ - ■ 

^~~~ z def«uJt_Frlnge :* 

M ax 4.000-03 @Elm 10557.1 
M in8. 500-05 <g>Elm 10327.1 

MSC.Patran 2005 r2 17- Oct-07 08:55.07 002639 

Fringe Derived Results. Subcase 4. Strain Tensor ^SCAPCL . . . (N ON-LAYERED) 

’• ’ 

^ m ~ 

t \ M JL1 

^ v 

| V 000227. 

000054 

default.Fringe 
Max .002639 @Elm 10659 1 



Min .000054 <$Eltn 10327.1 


• In original design the strain field displayed few small areas of high strain (exceeding 
allowable) with the rest of the field having relatively low strain values 

• Optimization process more evenly distributed the strain field 
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Run Time 


• CometBoards was run on a Red Hat Linux 2.6.9-67.ELsmp O/S, with 
x86_64 architecture, 2600 MHz, 4 cpus, 8GB of memory, 32 bit numeric format 

• MSC/Nastran version 2005.5.0 (2005R3) 

• FPI of Nessus level 6.2, dated 29 Sept. 1995 


Activity 

Time 

One analysis cycle - static 

5 sec. 

One analysis cycle - 20 frequencies 

51 sec. 

Deterministic optimization run (model- 1) 

39 min. 

Stochastic analysis (model -1 +load-A): 


- (303 1 number of random variables: 3025 ply strains + 6 displ) 

47 min. 

Stochastic optimization runs (GRC model + GRC load): 


- (61 number of p-levels) ~ (5 days + 8 hrs) 

128 hrs 

Stochastic optimization runs (LaRC model + GRC load): 


- (61 number of p-levels) ~ (5 days + 6 hrs) 

126 hrs 
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Summary 

• A single winner from the three methods can not be identified for design 
calculations 

- Traditional design is based on fully utilized concept. Weight is back calculated. 

- Deterministic optimization minimizes the weight with constraints on stresses and/or 
displacements - widely accepted method 

- Stochastic design similar to deterministic but parameters are random variables with 
specified reliability or 1 failure to N samples 


• All three methods produced acceptable solutions with some variation which 
may be due to the indeterminacy that can influence stress or strain flow in 
the members 

• Variation of weight was modest for all methods. In design variables up to 
one third variation was noticed 

• In stochastic design weight versus reliability traced out an inverted S -shape 
graph. Weight increased when risk was reduced and vice-versa 

• Stochastic design should be used with caution because the distribution 
function, mean value and standard deviation were assumed random 

variables 
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